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The structure of the equation of state lu could be very complicate in nature while a few linear 
models have been successful in cosmological predictions. Linear models are treated as leading 
approximation of a complete Taylor series in this paper. If the power series converges quickly, one 
can freely truncate the series order by order. Detailed convergent analysis on the choices of the 
expansion parameters is presented in this paper. The related power series for the energy density 
function, the Hubble parameter and related physical quantities of interest are also computed in this 
paper. 
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Recently there have been advances in our abilities in cosmological observations for the quest of exploring the 
expansion history of the universe. It carries cosmology well beyond determining the present dimcnsionlcss density of 
matter fl m and deceleration parameter go Q- It seems possible to reconstruct the entire function a(t) representing 
the expansion history of the entire universe. Earlier on, cosmologists sought only a local measurement of the first two 
derivatives of the scale factor a, evaluated at a single time to . Attempt will be made trying instead to map out the 
function determining the global dynamics of the universe in the near future. One notes that many qualitative elements 
\ of cosmology follow directly from the structure of the metricQ- Therefore, deeper understanding of our universe 

■ requires knowledge of the dynamics of the scale factor a(t) echoing the transition of energy between components from 
the epoch of radiation domination to that of matter domination. This is also known to be a key element in the growth 
of density perturbations into large structure. Indeed, many proposed observation will be able to probe the function 
a(t) more completely throughout all ages of the universe 0, 0, 0| - 

A number of promising methods are being developed including the magnitude-redshift relation of Type la super- 
novae. The goal of mapping out the recent expansion history of the universe is well motivated. The thermal history 
of the universe, extending back through structure formation, matter-radiation decoupling, radiation thermalization, 
primordial nucleosynthesis, etc. is very important in the study of cosmology and particle physics, high energy physics, 
neutrino physics, gravitational physics, nuclear physics, and so on 

Expansion history of the universe is similarly a very promising research focus with the discovery of the current 
^Jy acceleration of the expansion of the universe. This includes the study of the role of high energy field theories in 
the form of possible quintessence, scalar-tensor gravitation, higher dimension theories, brane worlds, etc in the very 
recent universe. The accelerated expansion is also important to the possible fate of the universe It is hence 

important to obtain the magnitude-redshift law relating to the scale factor-time behavior a(t) from these supernova 

■ observations with the proposed Supernova/ Acceleration Probe mission [Tl| . 

Therefore, the study of modelling different equation of state (EOS) derived from different theories plays an important 
role in the study of the recent expansion history of our universe. The study of dark matter and possible contribution 
from different combinations of different theories including the SUGRA and string theory will become important also in 
the near future once the LHC experiment starts. Therefore all practical methods to study its impact on the evolution 
history of our universe deserve full attentions. In fact, there have been proposals to parameterize different candidate 
models that are helpful in extracting important physics with the help of several different linear parameterizations of 
parameter z/(l + z) or z for the corresponding equation of state function u> 0, El EL El El El Hi IH3 ■ These 
linear models are shown to be very successful in recreating various physical functions. 

Linear models, e.g. lu = ujq + toiy with y = z/(z + 1), should be thought of as leading approximation of a complete 
Taylor series. If this is so, result derived from higher order in y, e.g. 0(y 2 ) terms, should only provide marginal 
contribution and hence can be ignored. If the higher order contribution is appreciable, one should be very careful in 
dealing with the truncated series. As a practical analysis, one should compute all related physical functions also as a 
power series order by order to obtain reliable observable according to the precision requirement. Otherwise, unphysical 
contributions may build up and leave the final result invalid. Fortunately, this complication may be unnecessary if 
the higher order contributions are small as compared to the leading order contribution. 

Therefore, we try to extend the parametrization further to include the effect of higher orders of the series expansion 
results. Related physical quantities are treated carefully order by order in order to extract more reliable information 
from these expansions. One also tries to determine the expansion coefficients from the fitting of the measured Hubble 
parameters or energy density. As a result, one may reconstruct the series expansion of the equation of state and probe 
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the nature and origin of the matter sources. 

Finally, we will also present an error analysis to find the range of convergence and possible error control for a 
meaningful truncation. For example, the linear model given by lu = —0.82 + 0.58y fits the SUGRA result to a good 
precision at large z even the linear model is already off 27% when z ~ 1.7^i|. One can show that the next leading 
expansion coefficient with us^y 1 < O.ly 2 only provides at most 20% error even when the redshift is close to 1100. As 
a result, one can show that a reliable truncation is possible when the higher order expansion coefficient is small for 
all range y = z/(z + 1) = — > 0.9991 up to the last scattering surface when z ~ 1100. Therefore, it does not matter 
much whether one sums up all power of y coupled to the leading coefficients luq and ux\ . In addition, the linear model 
to = —0.4 + O.lly approximates the inverse power law model very well at small y[l9|- Similar analysis will also be 
shown in this paper for comparison. 

In addition, the evolution of the equation of state lu of the physical universe could be very complicate. Indeed, 
the equation of state can be described by lu = p/p = (0 2 /2 — V)/(4> 2 /2 + V) for many different effective theory with 
effective potential V and scalar field <j>. For example, Vinv ~ 4>~ a an d VsuGRA ~ 0~ a exp[</> 2 ] for inverse power 
law models and SUGRA models respectively^^. The nature of our physical universe may also consist of many other 
combinations of scalar fields contributions. It is therefore important for us to figure out which model plays the most 
important role in the expansion history of our universe. The linear model or leading approximation appears to be of 
help in determining the essential part of the nature of the equation of state and its corresponding effective theory. 
Therefore a more detailed analysis on the convergence properties of various approximated models is very important 
as a research topic. 

Since one is not sure about the nature of the function u, the Taylor series expanding uj(y) — J2 n ^nU 71 into series 
sum with expansion parameter u n may also be used to determine leading expansion parameter with the help of future 
measurement. By fitting the Hubble parameter measurements one could determine leading terms of the expansion 
coefficients w n . As a result, one would be able to reconstruct the recent evolution of the function lu to the better 
precision. We will present the complete analysis in this paper. Since future observations will be able to provide tight 
constraint for only a few parameters in the near future, the material provides in this paper can be used to determine 
these parameters in a more reliable details. Hopefully, future development of observation tools may provide deeper 
insight to determine more constraints in these theories. 

We will present, in section II, a Taylor series expansion of some physical models. One will explain the choices of the 
expansion variables which come from the integration involving dz/(l + z). The convergence range of these expansion 
variables will also be presented. In section III, we will also compute leading terms of the energy density function p, the 
Hubble parameter H, and the conformal time rj for these expansions in order to discuss the difference between these 
approaches. (cf. 01 an d Fig. 1 of [l9j)- Due to the fact that the parameter z/(l + z) has a larger naive convergent 
range, this parameter could be more useful in large z expansion. We also compute and analyze the convergence range 
of a linear model in section IV. 



II. TAYLOR SERIES OF THE EQUATION OF STATE 

The equation of state is defined by the relation lu = p/p. In addition, the field equations of the universe can be 
shown to be: 

dlnp(z) = 3(w + l)dln(l + z), (1) 
& = *£p, (2) 

for a flat FRW universe with 1 + z the redshift defined via a/ao = 1/(1 + z). Various linear models [l^, have 
been suggested for the EOS in the literatures. We will try to point out in this paper the relations between these linear 
models and go further to obtain a more complete leading-order expansions. 

Taylor expansion is known to be one of the best ways to extract leading contributions from a generic theory with 
the help of a power series expansion of some suitable field variables. The Taylor series is normally convergent quickly 
depending on the structure of the expansion coefficients. The power series will, however, converge quickly if the range 
of variable is properly chosen. Hence it is important to choose an appropriate expansion variable for the purpose of 
our study. 

For example, one may expand the EOS, assuming to be a smooth function for all z, as a power series of the variable 
z around the point z = 0. This will lead to the power series expansion: 

n— 
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Here the summation of n runs from to oo and f( n '(z) — d n f(z)/dz n for any function f(z). From now on, the range 
of summation will be omitted for convenience throughout this paper unless it is different from the range from to oo. 

One sees clearly that there are three useful and practical ways to expand the equation of state: z, 1/(1 + z), and 
ln(l + z), from the structure of the conservation law given by Eq. (Q. It turns out that the expansion parameters 
will become z, —z/(l + z), and ln(l + z) respectively if one tends to expand the function 10 about the point z = 0. In 
practical, we will expand the second choice with the parameter z/(l + z) instead of —z/(l + z) for convenience. 

In this section, one will first expand the EOS as a power series of the variable w = 1/(1 + z) = a/ao around the 
point w(z — 0) = 1. The result is 

- £ ^^(, - ir - v ^^(-)"(t^)". m 

n n 

Note that if we expand the EOS as a power series of the variable y — z/(l + z) around y(z = 0) = 0, we will end 
up with a similar power series: 

u(y(z)) = 2^ - { y ■ (5) 

n 

The Taylor series are normally convergent quickly depending on the structure of the expansion coefficients. Nonethe- 
less, one would prefer to choose a more appropriate expansion parameter in order to make the series converge more 
rapidly. As a result, leading order terms will be enough to extract the most important physics from the underlying 
theory. Therefore the advantage of the y expansion is that the power series converges rapidly for all range of the 
parameter —1 < y — z/{\ + z) < 1, or equivalcntly, —1/2 < z < oo as compared to the range \z\ < 1 for the power 
series expansion of z shown in Eq. 0. 

One notes, however, that one should also expand all physical quantities and the field equations to the same order 
of precision we adopted for the EOS expansion. Higher order contributions will not be reliable unless one can show 
that the higher order terms does not affect the physics very much. For example, the liner order is good enough for the 
expansion of the EOS modelling SUGRA model [Rj. This is because that the linear term fits the predicted EOS for 
SUGRA model to a very high precision. One readily realizes, however, that it is not easy to track the series expansion 
order by order due to the form of the Eq. for the EOS for the y expansion. This is because that after performing 
the integration, one needs to pay attention to the distorted integration result. Indeed, one needs to write y — 1 — w in 
order to perform the integration involving d\nw. Indeed, one will need to recombine the result back to a power series 
of y. The trouble is that the lower order terms could sometimes hide in the higher order terms in w. It may not be 
easy to track clearly the lower order y expansion when we perform the expansion, for example, due to the exponential 
factor in the Eq. QJ. Therefore, one finds that the most natural way to expand the EOS is to expand it in terms of 
the variable x = ln(l + z) around the point x(z = 0) = 0. Indeed, this power series can be shown to be: 

«(*(*)) =£^*». (6) 

/ — 1 TV. 

n 

And this power series converge very rapidly in the range —1 < ln(l + z) < 1, or equivalently, —0.63 ~ 1/e — 1 < z < 
e — 1 ~ 1.72. Here e ~ 2.72 is the natural factor. Note that this limit happens to agree with the proposed scope of the 
SNAP mission. We will study these different power series expansion for the EOS and its applications in the following 
sections. 



III. SOME PRACTICAL EXPANSIONS 



Due to the structure of the differential equations Q , one finds that it will be easy for us to compute related Taylor 
series by expanding to either as functions of z, or y = z/(l + z), or x = In z. Therefore, we will present details of these 
expansion series in this section. For convenience, we will use repeated notation for the expansion series in this section 
for convenience and economics of notations. One should bear in mind that these coefficients are defined differently 
associated with different arguments y, x and z defined in each subsection. 



A. Power series of z/(l + z) 



Note that the expansion parameter y = z/(l + z) = 1 — a/a a . Therefore, the Taylor series expansion around small 
y is equivalent to a series expansion around a — ap. This is again a series expansion around the very recent universe 
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near a — clq. The expansion series in y has a naive convergence range for all y < 1 which is defined to be the early 
universe with a < do- We will show in details how to extract the leading terms in the y-expansion with y = z/(l + z) 
around the point y(z — 0) = 0. One can write the expansion coefficient as u) n = (1 + uip n '(y = 0)/n! such that the 
power series for the expansion of the EOS becomes 



l+uj(y(z)) = ^2u n 



y 



(7) 



Hence the Eq. can be shown to be 



dlnp(y) = 3(1 + u) 



dy 
1-1/ 



J2 



n+k+l 



n + k + 1 



(8) 



Note that we are now expanding with respect to the smooth function (1 + u>), instead of u), as a Taylor series for 
convenience. Therefore, one can integrate above equation to obtain 



P(y) = Po exp[ 3 + k + 1 ^ " P ° X ^ = Po 12 X " 



n+k+l 



(9) 



Note that one needs to expand the function p as a power series of y too in order to extract the approximated solution 
with appropriate order. The expansion coefficient X n is defined as X n — X^ n \y — 0)/n!. Here, the superscript in 
X 1 denotes differentiation with respect to the argument y of the function X(y). One can show that X' = XY with 
Y = 3 J2 n Sfc UJ n y n+k - In addition, one can show that 



(10) 



Hence one has 



Y«Hy = o) = 3(U) Y^n- 



(11) 



n=0 



One can also show, for example, that 



X" = X{Y 2 +Y') 
X'" = X(Y 3 + 3YY' +Y") 
X {i) = X(Y 4 + 6Y 2 Y' + 3(Y') 2 + AYY" +Y' 1 



(12) 
(13) 
(14) 



This series does not appear to have a more compact close form for the multiple differentiation with respect to y. One 
can, however, put the equations as a more compact format: 



X^+V =X[Y+— } l Y. 

dy 



(15) 



It appears, however, that one needs to do it manually even it is straightforward. We will only list the leading terms 
as this is already suitable expansion for our purpose at this moment. 
Hence one has 



*o = 1, 

Xi = 3cj , 

X 2 = -[9wg + 3o;o + 3wi], 

X 3 = -[9wo + 9w (wo+wi) + 2(wo + wi+c<^)], 

X 4 = -[27ujq + 54wo(wo + w\) + 9(ui n + uji) 2 + 24cj (^o + £*>i + u 2 ) + 6(o;o + uj\ + uj 2 + W3) 



(16) 
(17) 

(18) 
(19) 
(20) 
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Therefore, one can expand the final expression for the energy density p accordingly. Indeed, the result is 

p = p {l + 3uj y + -[9wq + + 3wi]y 2 + - [9wq + 9cj (^o + wi) + 2(w + wi + u 2 ))y 3 

+ 1[27uj% + 54u^(w + wi) + 9(w + ujx) 2 + 24w (cj + wi + u 2 ) +6(w + wi + w 2 + w 3 )]y 4 } + 0(y 5 ) 
8 

3 1 

= po{exp[3w y] + 2 I w o + wi]y 2 + -[9uj q (ujq + wi) + 2(cj + wi + w 2 )]y 

+ i[54u 2 (cj + wi) + 9(w + ujx) 2 + 24w (w + wi + uj 2 ) + 6(w + uji + uj 2 + uj 3 )]y 4 } + 0{y 5 ) (21) 
8 

to the order of y 4 . Note that one keep the order of precision to y 4 in computing the energy density p even we are 
expanding the EOS only to the order of y 3 . This is due to the special structure in the energy momentum conservation 
law 0). In addition, the largest power terms of ujq in the series come from the combination X)nP^"(0)/ n ']y™ = 
exp[3(jJo?/] can be summed over directly. Moreover, one can also show that the Hubble parameter H = Ht^X 1 ^ 2 with 
Ho = y/SnGpa/'i. And the expansion for X 1 / 2 can be obtained by replacing all u> n with uj n /2 in writing the expansion 
for X. Therefore one has 

3 3 1 1 

H = H {exp[-u> y] + -(u) Q +uji)y 2 + -[9cj (cj + Ui) + 4(w +u\ + uj 2 )]y 3 + — [ 

27uj%(u + uji) + 9(oj + ujx) 2 + 24uj (uj +lo 1 + lo 2 ) + 12(w + lo 1 + uj 2 + uj 3 )]y 4 } + 0(y 5 ). (22) 
Note also that one can also compute the conformal time according to the expression: 

which comes from the definition dq = dt/a. Knowing that 1/(1 — y) 2 = J2 n ( n + !)£/"" , one can show that 



HoV= I" dy'X-? V(n + 1) 
Jo 



y ■ (24) 



One can write X x l 2 = J2 k xuy k = J2 k Xk(u) n — > ~^n/2)y k , and expand H n — J2k 7 lky k for convenience. Therefore, 
one has 

En + 1 
— j — Zj-n-l (25) 

for Z > 1. Note that 770 =0. As a result, one can easily reconstruct the power series of Hot). Therefore, one has 

Hon = y - [^uj - l]y 2 + [^w 2 - - ^uj + l]y 3 - [-^uj 3 - ^uj 2 - J^o^i + + ^1 + - % 4 + 0{y 5 ) 
V 3 2 , r 3 2 1 5 3 9 3 27 2 9 13 1 ! 



—LUnV + \ — LU n Wi OJo\V — \ Wn W„ WnWi H OJn H OJi H LUr>]v + 0(ll ). (26) 

l-y 4 )y L 8 4 4 0JJ L 64 32 32 8 2 8 " J,y VW y V ; 

Note that terms independent of ojq, y + y 2 + ■ ■ ■ , comes from the combination K$ = YliLx VlV 1 = Eti 1 ( n + 
X)xi- n -iy /I. Taking the terms with n = I — 1 in n-summation, one can show that this summation becomes 

Ei=i ^o?/' = E/=i y l = y/0- - v)- 



B. power series of ln(l + z) 



Note that the expansion parameter x = ln(l + z) = — ln(a/ao). Therefore, the Taylor series expansion around 
small | a; | is equivalent to a scries around a = do- This is again a series expansion around the very recent universe near 
a = clq. Indeed, the expansion series in x has a naive convergence range for all < \x\ < 1 which is defined to be the 
early universe with a®/ e < a < a^. Note that we are interested only in the past universe where a < ao- We will show 
in details how to extract the leading terms in the x-expansion with x — ln(l + z) around the point x(z — 0) = 0. One 
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can write the expansion coefficient as u>„ = (1 + u)( n '(x — 0)/n! such that the power series for the expansion of the 
EOS becomes 



1 + Uj(x(z)) = Y^ L) n r j 



(27) 



Note that we use the same notation for u> n in different parameterizations for convenience. Hence the Eq. Q can be 
shown to be 



d In p(x) = 3(1 + ui)dx = d 



3w„ — — 
n + 1 



(28) 



Note that we are now expanding the physical quantities with respect to the function (1+lu) instead of w for convenience. 
Therefore, one can integrate above equation to obtain 



ii+i 



p(x) = pa exp[ 3 ] = p X(x) = p Y] ^n«™ 

* — ' n + 1 * — ' 



(29) 



Note that one needs to expand the function p as a power series of x too in order to extract the approximated solution 
with appropriate order. The expansion coefficient X n is defined as X n = X^ n \x — 0)/n\. One can show that 
X' = XY with Y = 3 J2n = 3 (! + w )- Therefore, one has 



One can also show, for example, that 

X" = 
X'" = 



In addition, one can show that 



X { 



y(0( x = o) = 3(11) uji. 



X(Y 2 +Y') 

X(Y 3 + 3YY' + Y") 



X {i) = X(Y 4 +6Y 2 Y' + 3(Y') 2 + 4YY" + Y'' 



n k 



(n + k-l + l)l 



(30) 



(31) 
(32) 
(33) 



(34) 



Therefore, one has 



i-i 



xW(0)=Y^3X„o; i _ n _ 1 (;-l)! 



(35) 



Hence one obtains the recurrence relation for the expansion coefficients of X ri 



i-i 



Xi — y X n Wi— n —\. 



(36) 



As a result, one has, for example, 



X = 1, 
Xi = 3w , 



X 2 
X 3 



1 



[Qug + awi], 

[9wq + 9tt> ^i + 2lu 2 ] 



(37) 
(38) 

(39) 
(40) 
(41) 
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Therefore, one can expand the final expression for the energy density p accordingly. Indeed, one has 
p = po{l + Sujqx + - [9ojq + 3u!i]x 2 + - [9^q + 9ujquji + 2cj 2 ]a; 3 

+ I [27w 4 + 54w^i + 9uo\ + 24w w 2 + 6w 3 ]a; 4 } + 0(x 5 ) 
8 

= po \ exp[3uj x} + \uix 2 + i[9w ^i + 2uj 2 ]x 3 + i[54w 2 u>i + 9uj\ + 24w w 2 + 6w 3 ]x 4 l + 0(x 5 ). (42) 



In addition, one can show that the Hubble parameter H = HqX 1 / 2 with H — ^8nGpo/3. And the expansion for 
X 1 / 2 can be obtained by replacing all to n with uj n /2 in writing the expansion for X. The result is 

f 3 3 1 1 1 

H = H < exp[-w a;] + -j^ix 2 + -[9w ^i + 4u> 2 ]x 3 + — [27w 2 Wi + 9oj\ + 24w w 2 + 12w 3 ]x 4 \ + 0{x 5 ). (43) 
1^ 2 4 8 32 J 

Note also that one can also compute the conformal time according to the expression: 

H Q r} = dz'X-* = dx'X-^- expfx'l = V / dx' X~^ — . (44) 
Jo Jo „Jo n - 

Therefore, one can easily compute the expansion of rj in a straightforward manner. One can write X^ 1 ! 2 = J^ fe XkX k = 
J2k A"fc(w n — * — u n /2)x k , and Hqij = J^k VkX k for convenience. Therefore, one has 



n=0 

for all I > 1. Note that r)o — 0. As a result, one can easily reconstruct the power series of HqT]. Indeed, one has 

Hoi] = x + i(2 - 3u )x 2 + 7^(9^o ~ 6cj i _ 12lj o + 4)x 3 

- Jg^( 27w o ~ 54w 2 - 54u wi + 36w + 36ui + 24w 2 - 8)x 4 + 0(x 5 ) 

3 1 1 

= e x - 1 - -w x 2 + -(3w 2 - 2cji - 4w )x 3 - — {9uj% - 18w 2 - 18w wi + 12^ + 12w a + 8cj 2 ).t 4 + 0(x 5 X-46) 

4 8 64 

Note that a trivial complete sum of terms x + x 2 /2\ + x 3 /3l ■ ■ ■ = e x — 1, derived from the n = 1 — 1 terms in Eq. H45JI. 
has been summed over for convenience in the final expression. 



C. power series of z 



Note that the expansion parameter z = (ao/a) — 1. Hence the range of z of interest is the range z < 1 in the past 
universe which is equivalent to the range where ao/2 < a < a . Therefore, the Taylor series expansion around small 
\z\ is equivalent to a series around a = ciq. This is again a series expansion around the very recent universe near 
a = clq. We will show in details how to extract the leading terms in the ^-expansion around the point z = 0. One can 
write the expansion coefficient as ui n = (1 + U)p n '(z = 0)/n! such that the power series for the expansion of the EOS 
becomes 

l + w(z) = Y,Un* n - ( 47 ) 

n 

Note that we use the same notation for u> n in different parameterizations for convenience. Hence the Eq. Q can be 
shown to be 

(48) 



dlnp(z) = 3(l+a;) 



dz 
l + z 



£(-l)*3w r 



z 



n+k+l 



n + k + 1 
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Note that we are now expanding the physical quantities with respect to the function (1+w) instead of u> for convenience. 
Therefore, one can integrate above equation to obtain 



p(z) = po exp 



3 y ( _l)^ 



y n-\-k-\-l 



n + k + l 



Po 



X(z) = P0 J2 X nZ n - 



(49) 



Note that one needs to expand the function p as a power series of z too in order to extract the approximated solution 
with appropriate order. The expansion coefficient X n is defined as X n — = 0)/n!. One can show that 

X' = XY with Y = 3J2„ k (-l) k uJ n z n+k . Therefore, one has 



y(0 ( 



3 l^^- 1 ) {n + k -l)^ Z 



n k 



Hence one has 



y«(z = o) = 3 i\ 



n+l 



n=0 



One can also show, for example, that 



X" = X(Y 2 +Y') 
X'" = X(Y 3 + 3YY' + Y") 
X {i) = X{Y i + 6Y 2 Y' + 3(Y') 2 + AYY" +Y r 



(50) 



(51) 



(52) 
(53) 
(54) 



This series does not appear to have a more compact close form for the multiple differentiation with respect to z. One 
can, however, put the equations as a more compact format: 



= X[Y + ±] l Y. 
dz 



(55) 



It appears, however, that one needs to do it manually even it is straightforward. We will only list the leading terms 
as this is already suitable expansion for our purpose at this moment. 
Hence one has 



X 
Xi 

x 2 

X; 
Xi 



1, 

3wo 
1 

2 
1 

2 
1 



[9lo — 3loq + 3lji], 

[9uJq + 9w (wi - ujq) + 2(uj — tJi + w 2 )], 



: [27wq + 54wq(^i - wo) + 9(w - uji) 2 + 24w (w - uj\ + lu 2 ) + 6(w 3 - lj 2 + Wi — ujq)]- 

a 

Therefore, one can expand the final expression for the energy density p accordingly. Indeed, the result is 
P 



(56) 
(57) 

(58) 
(59) 
(60) 



3 1 

p {exp[3w z] - -[oj - uj x )z + - [9w (wi - w ) + 2(uj - u\ + uj 2 )]z 3 



1 



+ -[54cj (tJi - lu ) + 9(luq — uji) + 2Alu (uj — ui + u 2 ) + 6(w 3 - lu 2 + oji — uj )]z } + 0(z ) 



(61) 



to the order of z 4 . Note that one keep the order of precision to z 4 in computing the energy density p even we are 
expanding the EOS only to the order of z 3 . This is due to the special structure in the energy momentum conservation 
law CJ. In addition, one can show that the Hubble parameter H = HqX 1 ^ 2 with H = ^SnGpo/S. And the 
expansion for X 1 / 2 can be obtained by replacing all uj n with w ra /2 in writing the expansion for X. Therefore one has 

3 3 1 

H = H {exp[-uj z} + -(wi - lu q )z 2 + -[9w (^i - ^o) + 4(w - ^i + uj 2 )]z 3 



+ — [27uj 2 (lu 1 - (j ) + 9(lu - loi) 2 + 24w (w - wi + w 2 ) + 12(cj 3 - w 2 + wi - lu )]z 4 } + 0(z 5 ). (62) 
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Note also that one can also compute the conformal time according to the expression: 

H r]= f dzfX-*. (63) 
Jo 

One can write X -1 / 2 = J2 k XkZ k = ^2 k Xk(w n — > —Lu n /2)z k for convenience. Therefore, one has 

oo 

^ = E^ fc ( 64 ) 

n=l 

As a result, one can easily reconstruct the power series of H rj. Therefore, one has 

H V = z- '^-z 2 + l[3uil + 2uj - 2wi]z 3 - ^t[9w^ + 18w 2 - 18w wi + 8(w - Wi + w 2 )]z 4 + 0{z b ) 
4 8 o4 

= g^(l - cxp[-^z] ) + I[ Wo - Wl ]z 3 - ^[3w 2 - 3w W! + 4(w - wi + w 2 )]z 4 + <3(z 5 ). (65) 

IV. COMPARISON AND ANALYSIS 

In summary, one has obtained the scries sum of the corresponding energy density, Hubble parameter and conformal 
time for different parameterizations of the function ui. For y = z/(l + z) one has: 

3 1 

P{y) = /°o{exp[3a> J/] + ^[^o + ^i]y 2 + -[9w (w + Wi) + 2(w + uni + w 2 )]y 3 

+ -[54wg(w + wi) + 9(w + ^i) 2 + 24w (w + wi + w 2 ) + 6(w + wi + w 2 + w 3 )]y 4 } + 0(y 5 ) (66) 
8 

3 3 1 1 

H{y) = iJo{exp[-w y] + -(w + ^\)y 2 + -[9w (w + ^i) + 4(w + wi + w 2 )]y 3 + ^[ 

27w 2 (w + wi) + 9(w + wi) 2 + 24w (w + wi + w 2 ) + 12(w + <*>i + ^2 + w 3 )]y 4 } + 0(y 5 ). (67) 
3 3 15 9 27 9 13 1 1 

H v(y) = - ^oy 2 + [g^o - - | w o]y 3 - - 3^0 - 32^1 + y^o + + g^ 2 ]y 4 + O(y 5 0S8) 

For x = ln(l + z), one has: 

p(x) = po jexp[3w x] + ^wix 2 + i[9w ^i + 2w 2 ]a; 3 + ^[54w 2 wi + 9w 2 + 24w w 2 + 6w 3 ]a; 4 j + 0{x b ). (69) 

{331 1 1 

cxp[-w £] + + ^[9w ^i + 4w 2 ]a; 3 + — [27w 2 wi + 9w 2 + 24w w 2 + 12w 3 ]a; 4 \ + 0(x 5 ). (70) 

2 4 8 32 J 

3 1 1 

H r](x) = c x - 1 - -w x 2 + -(3w 2 - 2wi - 4w )x 3 - — (9w 3 - 18w 2 - 18w wi + 12w + 12wi + 8w 2 )x 4 + 0(a;f|l) 

4 8 64 

For the 2 expansion, one has: 

3 1 

p(z) = po{cxp[3cj z} - -[uj - uji]z 2 + -[9w (wi - w ) + 2(w - u>i + w 2 )]z 3 

+ i[54a) 2 (wx - w ) + 9(w -^i) 2 + 24w (w - wi + w 2 ) + 6(w 3 - w 2 +u 1 - u a )]z 4 } + 0{z 5 ) (72) 
8 

3 3 1 

H{z) = iJ {exp[-w z] + -(u>i ~w )z 2 + -[9w (wi - w ) + 4(w - Wi + w 2 )]z 3 
2 4 o 

+ -'-[27w 2 (w 1 - w ) + 9(w - ^i) 2 + 24w (w - wi + w 2 ) + 12(w 3 - w 2 + Wi - w )]z 4 } + 0(z 5 ). (73) 

flo»?(«) = ^(l-cxp[-^z] ) + i[w -w 1 ]z 3 -^[3w 2 -3w w 1 +4(w -w 1 +w 2 )]z 4 + O(z 5 ). (74) 

Note that the expansion coefficients u> n are defined differently for different expansions even we are using the same 
notation for convenience. Also note that at small red-shift z <C 1, y and x are both very close to z. This is the reason 
why we are ending up with similar leading term in each expansion. Once the red-shift z is extended, difference in the 
leading terms will be significant. One is therefore able to distinguish the contribution of the next-leading terms. 
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As mentioned earlier, one should be able to determine the expanding coefficients u n with the result from the 
measurements of future experiments. Once the expansion coefficients are determined from fitting with the Hubble 
parameter or energy density, one will be able to reconstruct the function to to higher order and higher precision. The 
result should be helpful to determine the origin and the nature of the matter sources and to be compared with some 
possible combinations of various fundamental theories. 

The linear parametrization of the to has been applied to the study of the evolution history of our universe. The 
linear models adopted has been shown to be useful in making predictions for future observations [l3| . One is naturally 
lead to answer the question whether one should take the linear model as a complete theory and integrate it without 
worrying about the higher order corrections. To be more specific, linear model like lu(z) = ui' +uo\z has been taken as 
a complete theory. Related physical quantities are computed with all power of u>\, in the power series shown above, 
summed over for the final result. On the other hand, if the linear model is taken as leading order approximation, one 
should ignore higher order contributions from z™(Vn > 2). 

If higher order corrections is small and negligible, it does not matter much whether one ignores or includes all higher 
order (in u>i) corrections. For example, in a convergent series, higher order contributions are smaller and smaller order 
by order. Therefore, it is fine to sum over all terms related to the expansion parameter u)\. Otherwise, one should 
pay attention to the convergent properties of these expansions for possible deviations. 

Looking at our results, it is easy to find that a re-summation of the expansion coefficients oj n in these formulae 
is difficult to obtain except the re-summation of loq which has been already shown in above equations. Fortunately, 
one can do it by a different way which will be shown in a moment. As a result, one would be able to see whether an 
expansion in power of to n is reliable or not. 

In next section, we will show how to sum up all terms order by order as series of expansion coefficients io n . Taking 
the linear model to = lo' q + LO\y as a complete theory and include all effect of iv' and io\ include more effect from the 
higher order y n contribution as compared to the leading order approximation approach. The inclusion of these effect 
may not be reliable unless they are small as compared to the leading terms. 

In addition, taking linear model as a complete theory, the effects from all higher order in ll>„>i are ignored. The 
truncation may need to appreciable error in evaluation of related physical observables unless these higher order terms 
are also small. We will also analyze the error due to the contributions from the higher order oj n >i term. 

Note again that linear model may serve as a good approximation theory in y expansion because the naive convergence 
range y < 1 covers the entire history of our universe from z = — > oo. If the linear model fits the original theory 
very well at small y and higher order contributions are not appreciable, linear model can be extended to remain valid 
even at larger y. The higher order effect related to L02 will be evaluated also for the hope that future experiments 
may provide better resolution to distinguish possible minor effect. 

In addition, one more advantage of the linear model written as u) = lo' q + io\y over the linear model in expansion of 
z can be easily seen from the relation: 

Lu{y) = oJq + ijjyy = ui' Q + ui- = lo' + uiiz(l — z + z 2 — ■ ■ ■ ) = oj(z) — lo\z 2 {\ — z + z 2 + • • • ). (75) 

_1_ ~\~ z 

Indeed, one can see that linear model oo(y) contains a few more higher order z n>1 effect than the linear model 
lo(z) = lo' + lo\z. The effect of the higher order terms depends, however, on the actual deviation from the underlying 
theory like SUGRA or inverse power law models. Evidences show that linear model in y works better than the linear 
model in z at large redshift for both reasons. 

V. ERROR ESTIMATION OF THE LINEAR MODEL 

Note that the expansion parameter y = z/(l +z) = 1 — a/ao. Therefore, the Taylor series expansion around small y 
is equivalent to a series expansion around a = ao. This is again a series expansion around the very recent universe near 
a = ao ■ The expansion series in y has a naive convergence range for all y < 1 which is defined to be the complete early 
universe with a < ao- We will show in details how to extract the leading terms in the y-expansion with y = z/(l + z) 
around the point y(z = 0) = 0. 

Indeed, one can write the expansion coefficient as uj n = (1 + oj) {n) {y = 0)/n! such that the power series for the 
expansion of the EOS becomes 1 + uj(y(z)) — to n y n . One can write w = 1 — y and write it as 

oc 00 n 

l + u(y) = 5>„2/ n = co + £> n (l - w) n = w(0) + J2 J2 C k UJ n(~w) k . (76) 

n n—1 n—1 k—1 
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0.1 0.2 0.3 0.4 0.5 

FIG. 1: Thin dotted line plots pi(y)/po and thick solid lines represents p 1 (y)/po for the linear model ui = —0.82 + 0.58y. 



Here u(0) = Y^Lo^n = [1 + u {y)]y=o- Hence the Eq. Q can be shown to be 

dy 



dhxp{y) = 3(l + o;) 



. dw 

-3(1 + w) — 

w 



^3d 



n=lk=l 



n uj n {-wf 



Therefore, one can show that 



p(y) = p(0) 



oo n k 



(1 - y )-M"> eM-3j2J2H 



k^ny 



(77) 



(78) 



n=l k=l 1=1 

In fact, one can show that when y — > 1, the energy density shown in Eq. (|78|) is dominated and proportional to 

p(y -> 1) K (1 - y)- 3 ^ ) (79) 

implying p oc (1 — ^-3(o.i8+u>2+w 3 +"0 r or moc j e i approximated by uj = —0.82 + 0.58y. This term will diverge at 
y = 1 if all higher order coefficients are small as compared to the leading coefficients such that the sum a>(0) remain 
positive. Therefore, one can expand the density as a power series of to n from above equation. For comparison, one 
can show that 



P2(y) = Pi(y) 



y 

(1 - y)exp[y+ —] 



2 i —ZuJ2 



(80) 



keeping terms to the order of u>2- Here p\(y) denotes p(y) with uj n = for all n > 2. Similarly, p2(y) denotes p(y) 
with bj n = for all n > 3. Hence the last factor of above equation provides the u>2 corrections of the power series. 

Note that linear model lu = —0.82 + 0.58y is considered as a complete model with the corresponding energy density 
p = p\ shown in Eq. (|SU|) . As a complete model, all powers of y n are included in the evaluation of energy density p\. 
Indeed, one can show that 



Pi = Po(l - y)- 3 ^ 0+ ^ exp[-3u iy ], 

P 1 = Po{l + 3^o2/ + 2 1 9w o + 3w o + 3^i]y 2 } 



(81) 
(82) 



with p 1 given above denoting the truncated energy density to the order of y 2 for the approximated counterpart. Results 
shown in Fig. 1 indicates that the difference between p\ and p 1 is small when y is small. But it becomes appreciable 
at larger y. Note again that u>o = 1 + uj(y = 0) = 0.18 in our notation for the linear model w — —0.82 + 0.58y. 
Similarly, one can show that H 2 — 8irGp/3 can be integrated to give: 



H(y) - H(0) 



(l-,)^(o)/ 2exph 3^^^ c^^f) 



n=l k=l 1=1 



(83) 
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FIG. 3: The equation of state u> is plotted as a function of y. Lower thick solid line denotes the contribution of linear model 
lu = —0.82 + 0.58y. Lower thin dotted line represents the addition of the y 2 correction with wiy 1 = O.lj/ 2 . Similarly, Upper 
thick solid line denotes the contribution of linear model to = —0.4 + O.lly. Upper thin dotted line represents the corresponding 
y 2 correction from L02y 2 = O.lj/ 2 . 



Therefore, one can also expand the density as a power series of uj n from above equation. Indeed, one can show that 

-3w 2 /2 



H 2 (y) = H 1 (y) 



y 

(1 - y)exp[y + — ] 



(84) 



Here H\(y) denotes H(y) with oj n = for all n > 2. Similarly, H2(y) denotes H(y) with u> n = for all n > 3. Hence 
the last factor of above equation provides the L02 corrections of the power series. 

In Fig. 2, the error p 2 / 'pi — 1 due to the presence of the ui 2 correction, as compared to the linear model, is plotted 
in 3D format as a function of y and u) 2 . It is easy to see that the error becomes appreciable when y > 7 and uj 2 > 0.4. 

In Fig. 3, the equation of state u> is plotted as a function of y. Lower thick solid line denotes the contribution of 
linear model oj = —0.82 + 0.58y 19]. This linear model fits SUGRA equation of state very well at small y. Lower 
thin dotted line represents the addition of the y 2 correction with us 2 y 2 = O.ly 2 . In addition, the upper thick solid line 
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denotes the contribution of linear model u> — —0.4 + 0.11?/|l9|. This model fits the Inverse power law model very well 
at small y. Upper thin dotted line represents the corresponding y 2 correction from Lo^y 2 = O.ly 2 . Note that it has 
been shown from Figure 1 of Reference 0] , linear models shown above fits the SUGRA and inverse power law model 
very well at small y as compared to the linear model in z. u> as a 

It is easily seen that a small y 2 correction with L02 = 0.1 does not modify the function uj very much. In addition, 
the correction to the density function for 0J2 = 0.1 is not appreciable according to Fig. 2. Therefore, the linear model 
can be truncated at the order of u>i for small y without affecting the final result in an appreciable way. 

From above presentations, one also finds that there are two comparisons required for the reliability for the linear 
model. First of all, one should check if the contributions from higher order coefficients w„>i will affect the evaluation 
results. Secondly, one should also check if the inclusion of higher order terms y n>2 in the physical functions, e.g. pi, 
will affect the final result in an appreciable way. 

Due to the fact that the naive convergence range of y = z/(l + z) expansion covers the entire history of our 
early universe, this expansion appears to be the best way to approximate our physical models. As a result, the 
approximated linear model may represent the complete theory effectively for a larger domain y provided that the 
higher order corrections derived from 102 is small. Our results indicates that this appears to be the case for the 
SUGRA models and inverse power law models |l9|. Therefore, the corresponding linear model appears to be more 
reliable than the other linear models expanded in z or ln(l + z). 

Detailed analysis indicates that the linear model is a very reliable approximation within small y. Beyond the 
region of small y, one should be very careful with the higher order contributions when comparing with the future 
experiments. In practice, the sources of equation of state ui may be very complicate. In spite of the fact that it may 
be derived from complicate combinations of many different sources, linear model provides an easy way to determine 
the leading coefficients that should be very reliable at small y corresponding to our very recent universe. If one is able 
to determine the next leading term coefficients u>2 from better precision measurements in the near future, it would 
provide better way to distinguish the nature of the source of equation of state. 

VI. CONCLUSION 

Linear models, e.g. uj = cj' a + u)\y, should be thought of as leading approximation of a complete Taylor series. If 
this is so, result derived from higher order in y, e.g. 0(y 2 ) terms, should only provide marginal contribution and 
hence can be ignored. If the higher order contribution is appreciable, one should be very careful in dealing with the 
truncated series. As a practical analysis, one should compute all related physical functions also as a power series order 
by order to obtain reliable observable according to the precision requirement. Otherwise, unphysical contributions 
may build up and leave the final result invalid. Fortunately, this complication may be unnecessary if the higher order 
contributions are small as compared to the leading order contribution. 

Indeed, linear models of equation of state are considered as effective models for some underlying theoretical models. 
Many theories including Sugra and Inverse power law model can approximated by linear model described by lo = 
lu'q + u>xy. Therefore these linear models are potentially good candidate to simulate the large redshift properties 
of the underlying theories because the higher order corrections could be negligible at larger z, corresponding to 
y = z/(l + z) -> 1. 

Since z — > 00 is equivalent to y — ► 1. Hence y-expansion with y — > 1 is still in the range of naive convergence of 
the corresponding Taylor series. Therefore, as long as the higher order terms are negligible as compared to the liner 
approximation, one is free to consider the linear model as a well-behaved representation of the underlying theory. 
Otherwise, one should pay attention to possible deviations derived from the higher order corrections. 

We have calculated two possible errors when one treats the linear model as a complete theory instead of the leading 
approximation for an underlying theory. One possible error comes from the inclusion of all higher order terms (in y n 
) related to loq and io\ in the linear model calculation. The other error comes from the truncated higher order terms 
in u) n for the underlying theory. They are both computed and compared in section V. 

We have also tried to analyze possible deviations for the Sugra and inverse power law models in this section. Results 
show that linear model lo = lo' q + uj\y for these two models appears to be a good approximation even at large redshift 
region as long as the higher order corrections are small. 

Since we are only able to determine only a few leading coefficients in the future experiments, linear model appears 
to serve this purpose very well. Hopefully, one should, however, be able to determine for example the deviation due to 
the u>2 contributions and further the understanding of the underlying theory if better resolution can be made possible 
in the future experiments . 

In addition, we have also tried to provide a more complete list of the Taylor series in this paper. Hope that these 
presentations can provide better information for the quest of the mapping of our early universe. 
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Therefore, we try to extend the parametrization further to include the effect of higher orders of the series expansion 
results. Related physical quantities are treated carefully order by order in order to extract more reliable information 
from these expansions. One also tries to determine the expansion coefficients from the fitting of the measured Hubble 
parameters or energy density. As a result, one may reconstruct the series expansion of the equation of state and probe 
the nature and origin of the matter sources. 

As mentioned above, we also present an error analysis to find the range of convergence and possible error control 
for a meaningful truncation. The linear model given by to = — 0.82 + 0. 58y is known to fit the SUGRA result to a good 
precision at large z even the linear model is already off 27% when z ~ 1.7. One shows that the next leading expansion 
coefficient with uj2y 2 < O.ly 2 only provides at most 20% error to the density function even when the redshift is close 
to 1100. The error is even down to 10% for the Hubble function H. As a result, one show that a reliable truncation 
is possible when the higher order expansion coefficient is small for all range y = z/(z + 1) = 0^ 0.9991 up to the last 
scattering surface when z ~ 1100. Therefore, it does not matter much whether one sums up all power of y coupled to 
the leading coefficients ojq and uj\ as long as the higher order corrections due to u>2 and w„>2 are negligible. 

In spite of the fact that equation of state u> may be derived from complicate combinations of many different 
sources, linear model provides an easy way to determine the leading coefficients that should be very reliable at small 
y corresponding to our very recent universe. If one is able to determine the next leading term coefficients 0J2 from 
better precision measurements in the near future, it would provide better way to distinguish the nature of the source 
of equation of state. 

Indeed, the equation of state can be described by uj = p/p = {<j) 2 /2 — V)/(4> 2 /2 + V) for many different effective 
theory with effective potential V and scalar field <j>. For example, Vinv ~ < P^ a an d Vsugra ~ 4>~ a exp[^> 2 ] for inverse 
power law models and SUGRA models respectively The nature of our physical universe may also consist of many 
other combinations of scalar fields contributions. It is therefore important for us to figure out which model plays the 
most important role in the expansion history of our universe. The linear model or leading approximation appears 
to be of help in determining the essential part of the nature of the equation of state and its corresponding effective 
theory. Therefore a more detailed analysis on the convergence properties of various approximated models is very 
important as a research topic. 

In summary, one has evaluated the series expansion for the energy density, Hubble constant and the conformal time 
for the y, x and z expansion of the corresponding equation of state function uj up to order four in previous section. 
The proposed Supernova/ Acceleration Probe (SNAP) will carry out observations aiming to determine the equations 
of state of the energy density, providing insights into the cosmological model, the nature of the accelerating dark 
energy, and potential clues to fundamental high energy physics theories and gravitation. As a result, we show all 
suitable ways to parameterizing the equation of state for application to study its effect on the expansion history of 
the recent universe. 

A detailed discussion on the choices of the expansion parameters for the Taylor series of the equation of states 
uj is presented in this paper accordingly. For example, the Taylor series of the EOS is expanded as power series of 
the variables y — z/(l + z), x = ln(l + z) and z respectively. Due to the fact that the naive convergence range of 
y = z/(l + z) expansion covers the entire history of our early universe, this expansion appears to be the best way 
to approximate our physical models. As a result, the approximated linear model may represent the complete theory 
effectively for a larger domain y provided that the higher order corrections derived from UJ2 is small. This appears to 
be the case for the SUGRA models and inverse power law models. Therefore, the corresponding linear model appears 
to be more reliable than the other linear models expanded in z or ln(l + z). 

We also show how to obtain the power series for the energy density function, the Hubble parameter and related 
physical quantities of interest. The method presented here will have significant application in the precision distance- 
redshift observations aimed to map out the recent expansion history of the universe, including the present acceleration 
and the transition to matter dominated deceleration. 

Since we can power-expand all smooth EOS into a convergent power series for a reasonable range of the expansion 
parameters, it is more practical for the future probe to determine the expansion coefficients uj n or equivalently the 
local derivatives of the EOS. One may need to use different expansion series depending on the convergent speed of 
the power series. For example, it appears that the leading order term in the y expansion is better as a nice result 
for the SUGRA prediction. This is because the leading term is close enough to the theoretical prediction at small 
y region[13|. Nonetheless, one expects that fitting for a few more leading terms in the Taylor series will be able to 
provide us better information about the nature of the function uj. In addition, the result shown here is independent 
of the choice of the time tg. The local measurement of the expansion coefficients can be extended to the comparison 
of the expansion coefficients at any time. 
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